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ABSTRACT 

We study the axisymmetric, time-dependent hydrodynamics of rotating flows 
that are under the influence of supermassive black hole gravity and radiation from 
an accretion disk surrounding the black hole. This work is an extension of the 
earlier work presented by Proga, where nonrotating flows were studied. Here, we 
consider effects of rotation, a position-dependent radiation temperature, density 
at large radii, and uniform X-ray background radiation. As in the non-rotating 
case, the rotating flow settles into a conflguration with two components (1) an 
equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the 
system along the pole. However, with rotation the flow does not always reach 
a steady state. In addition, rotation reduces the outflow coUimation and the 
outward flux of mass and kinetic energy. Moreover rotation increases the outward 
flux of the thermal energy and can lead to fragmentation and time-variability of 
the outflow. We also show that a position-dependent radiation temperature can 
signiflcantly change the flow solution. In particular, the inflow in the equatorial 
region can be replaced by a thermally driven outflow. Generally, as it have 
been discussed and shown in the past, we flnd that self-consistently determined 
preheating/ cooling from the quasar radiation can signiflcantly reduce the rate at 
which the central BH is fed with matter. However, our results emphasize also 
a little appreciated feature. Namely, quasar radiation drives a non-spherical, 
multi-temperature and very dynamic flow. These effects become dominant for 
luminosities in excess of 0.01 of the Eddington luminosity. 
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1. Introduction 

The key property of Active Galactic Nuclei (AGN) is that they emit enormous amount 
of electromagnetic radiation over a very broad energy range. The AGN central location in 
their host galaxies imply that AGN radiation can play a very important role in determining 
the ionization structure and dynamics of matter not only near the AGN but also on larger, 
galactic and even intergalactic scales (Ciotti & Ostriker, 1997, 2001, 2007; Silk & Rees 1998; 
King 2003; Murray et al. 2005; Sazonov et al. 2005; Springel et at. 2005; Hopkins et al. 
2005; Wang et al. 2006; Fabian et al. 2006; Thacker et al. 2006, and references therein). 

In the first paper of this series (Proga 2007, Paper I hereafter), we reported on results 
from our first phase of gas dynamics studies in AGN on sub- and parsec-scales. This is 
a complex problem as it involves many aspects of physics such as multi-dimensional fluid 
dynamics, radiative and magnetic processes. Therefore, our approach was to set up sim- 
ulations as simple as possible and to start with exploring the effects of the X-ray heating 
[important in the so-called preheated accretion, (e.g., Ostriker et al. 1976; Park & Ostriker 
2001, 2007)] and radiation pressure on gas that is gravitationally captured by a black hole 
(BH). We adopted the numerical methods developed by Proga et al. (2000, PSKOO here- 
after) for studying radiation driven disk winds in AGN. Generally, our simulations cover a 
relatively unexplored range of the distance from the central BH: we end where models of 
galaxies begin (e.g., Giotti & Ostriker 2007; Springel et at. 2005) and we begin where models 
of BH accretion end (e.g., Hawley & Balbus 2002; Ohsuga 2007). 

In Paper I, we presented results from axisymmetric time-dependent hydro dynamical 
(HD) calculations of gas flows. The flows were non-rotating and exposed to quasar radiation. 
We took into account X-ray heating and the radiation force due to electron scattering and 
spectral lines. To compute the radiation fleld, we considered an optically thick, geometrically 
thin, standard accretion disk as a source of UV photons and a spherical central object as a 
source of X-rays (a corona). The gas temperature, T and ionization state in the flow were 
calculated self-consistently from the photoionization and heating rate of the central object. 

We found that, for a 10^ M© black hole with an accretion luminosity of 0.6 of the 
Eddington luminosity, the flow settles into a steady state and has two components: (1) an 
equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along 
the disk rotational axis. The inflow is a realization of a Bondi-like accretion flow. The second 
component is an example of a non-radial accretion flow which becomes an outflow once it 
is pushed close to the rotational axis where the radiation pressure accelerate it outward. In 
some cases the outflow is heated by radiation so that it can be accelerated also by thermal 
expansion. Our main result was that the existence of the above two flow components is 
robust to the outer boundary conditions and the geometry and spectral energy distribution 
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of the radiation field. However, the flow properties are not robust. In particular, the outflow 
power and collimation is higher for the radiation dominated by the UV/disk emission than 
for the radiation dominated by the X-ray/central engine emission. Our most intriguing result 
was that a very narrow outflow driven by radiation pressure on lines can carry more energy 
and mass than a broad outflow driven by thermal expansion. 

Here, we report on results from simulations that are basically reruns of those presented 
in Paper 1, but with inclusion of gas rotation. Our goal is to assess how rotation changes 
the flow solution. In particular, we study how rotation changes the flow pattern, mass and 
energy fluxes, and temporal behavior. We also present results from a new set of simulations 
that illustrate how complex and dynamic the flow evolution can be even for relatively simple 
initial and boundary conditions. We describe our calculations in Section 2. We present our 
results in Section 3. The paper ends, in Section 4, with discussion and our conclusions. 



2. Method 

In this paper we extend the work presented in Paper I by relaxing some assumptions 
and simpliflcations. Our numerical HD calculations are in most respects as described by 
Paper I. Here we only describe the key elements of the calculations and list the changes we 
made. We refer a reader to Paper I and PSKOO for details. 

We consider an axisymmetric HD flow accreting onto a supermassive BH. The flow is 
non-spherical because it is irradiated by an accretion disk. The disk radiation flux, J-'disk 
is highest along the disk rotational axis and is gradually decreasing with increasing polar 
angle, 9: J^disk |cos^|. The flow is also irradiated by a corona. We account for some 
effects of photoionization. In particular, we calculate the gas temperature assuming that 
the gas is optically thin to its own cooling radiation. We include the following radiative 
processes: Compton heating/cooling. X-ray photoionization heating, and recombination, 
bremsstrahlung and line cooling. We also take into account some effects of photoionization 
on radiation pressure due to lines (line force). Namely, we compute the line force using 
a value of the photoionization parameter, ^ and the analytical formulae due to Stevens & 
Kallman (1991). This procedure is computationally efficient and gives good estimates for 
the number and opacity distribution of spectral lines for a given ^ without detail information 
about the ionization state (see Stevens & Kallman 1991). Additionally, we take into account 
the attenuation of the X-ray radiation by computing the X-ray optical depth in the radial 
direction. To be consistent with our gas heating rates where we include X-ray photoionization 
but not UV photoionization, we do not account for attenuation of the UV radiation. 
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We assume that the total accretion luminosity, L has two components: Ldisk = fdiskL 
due to the accretion disk and = f^L due to the corona. For simplicity, we assume that 
the disk emits only UV photons, whereas the corona emits only X-rays i.e., the system UV 
luminosity, Lyy = fvvL = Ldisk and the system X-ray luminosity, Lx — fxL — (in other 
words /uv = /disk and /x = /*). 

With the above simplifications, only the corona radiation is responsible for ionizing 
the flow to a very high ionization state. In our calculations, the corona contributes to the 
radiation force due to electron scattering but does not contribute to line driving We note 
that metal lines in the soft X-ray band may have an appreciable contribution to the total 
radiation force in some cases. On the other hand, the disk radiation contributes to the 
radiation force due to both electron and line scattering. 

We perform our calculations in spherical polar coordinates (r, 6, (f) assuming axial sym- 
metry about the rotational axis of the accretion disk [9 = 0°). 

Our computational domain is deflned to occupy the angular range 0" < ^ < 90° and 
the radial range r[ = 500 < r < = 2.5 x 10^ r^,, where r^, = 3rs is the inner 
radius of the disk around a Schwarzschild BH with a mass, Mbh and radius rs = 2GMbh/c^- 
The r — 6 domain is discretized into zones. Our numerical resolution in the r direction 
consists of 140 zones. We fix the zone size ratio, drk+i/drk = 1.04 (i.e., the zone spacing is 
increasing outward). Gridding in this manner ensures good spatial resolution close to the 
inner boundary, ri. In the 9 direction, our numerical resolution consists of 50 zones. The 
zone size ratio is always d9i/d9i+i — 1.0 (i.e., grid points are equally spaced). 

For the initial condition, in Paper I, we assumed spherical symmetry and set all HD 
variables to constant values everywhere in the computational domain. Here we will allow 
the gas to rotate. 



In simulations with rotation, we break spherical symmetry of the initial and boundary 
conditions by introducing a small, latitude-dependent angular momentum. Namely, for large 
radii we assume the specific angular momentum, depends on the polar angle, 9, as 



2.1. 



Gas Rotation 
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where / = 1 on the equator {9 — 90°) and monotonically decreases to zero at the poles (at 
^ = 0° and 180°) . The initial distribution of the rotational velocity is: 

f for r < 10^ r 

'''^^'■'^^ " I //sin^r for r > 10^ rj 

We express the specific angular momentum on the equator as 

lo = \/r^/6 c r*, (3) 

where is the "circularization radius" on the equator in units of for the Newtonian 
potential (i.e., GM/r^ = v'^/r at r = r^r*). We adopt two forms for the function f{9): 

/■ ^flN _ / ^ < and 9 > 180° - 9^ 

' ~ \ for 9o< 9 < 180° -9o, ^ ' 



and 



/2(^) = l-|cos^|. (5) 



We also follow Paper I in setting the boundary conditions at the outer radius, Tq: to 
represent steady conditions at the outer radial boundary. During the evolution of each model, 
we apply the constraints that in the last zone in the radial direction, all HD quantities, expect 
the radial velocity, are set to their initial values at all times. The radial velocity is allowed 
to float. In this paper we set t'(/,(ro, 9) using eq. 2 whereas in Paper I, v^{ro, 9) was set to 
zero consistently with the initial conditions. 



2.2. Mean Radiation Energy 

One of the simplifications made in Paper I was an assumption that the radiation tem- 
perature, Tr [or more generally the spectral energy distribution (SED)] does not change with 
the position in the flow and it was set to the temperature of the isotropic corona radiation, 
Tx. This assumption finds justification in X-ray observations which show that quasar radi- 
ation heats a low-density gas nearly uniformally, on scales comparable to the Bondi radius, 
up to an equilibrium Compton temperature of about 2 x 10^ K (Sazonov et al. 2005; Allen et 
al. 2006 and references therein). However, within our theoretical framework, the SED and 
7r can change with position. For example, we assume that the UV/disk radiation decreases 
with increasing 9 (see eq. 6 in Paper I). Thus even if the X-ray/central object radiation 
does not change with 9 (see eq. 7 in Paper I) the ratio between the X-ray and UV fluxes 
increases with increasing 9. Consequently, our model predicts that radiation is softer near 
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the poles than near the equator. Moreover, we take into account the attenuation of the 
X-ray radiation but not of the UV radiation. Thus, for a given 9 the ratio between the X-ray 
and UV fluxes can decrease if the X-ray optical depth is high. To test how important these 
effects are, we allow Tr to vary with position. 

We introduce position-dependence of Tr, in the following way. We start with the stan- 
dard expression for the net Compton heating rate: 

£ = ne^(A;TR-4A;T), (6) 

where T is the radiation flux, rig is the electron number density (other symbols have their 
usual meaning). 

As metioned above, we consider two sources of radiation: a disk that emits UV photons 
with energy between and 50 eV, and a corona that emits photons with energy above 50 eV. 
Each of this sources has its own mean photon energy, < e > or equivalently radiation tem- 
perature, (< e >= /cTr). The radiation temperature of these two sources can be computed 
from: 

/■50 eV /-SO eV 

^Tdisk = / huj^disk,i.diy/ / J^disk,i.du (7) 

Jo eV Jo eV 



and 



kTx = / hfTy^^ydv/ I J^x,udv- (8) 

^50 eV JbO eV 

We consider these temperatures as additional free parameters. Wc chose the radiation tem- 
perature of the disk to be 2 x 10^ K and of the corona to be 2.9 x 10^ K. 

Having set the radiation temperature of the disk and corona radiation we can compute 
the radiation temperature of the total radiation field from: 

rj. rj. fxRT exp(-rx) + 2/uv cos 9 
fx exp(-Tx) + 2/uv COS 9 

where Rt = Tx,max/2uv. To obtain the expression above we used the definition of Tr and 
the formulae for the radiation fiuxes from the disk and corona (eqs. 6 and 15 in Paper I). 
Practically, to account for the position-dependent radiation temperature, we adopt the same 
formulae for the radiative heating/cooling rates as in Paper I but replace Tx with Tr (see 
eqs. 19 and 20 in PSK). The seventh column in Table 1 shows the values of the adopted Tr 
or its range for cases when we use eq. 9. 

We finish this section with a note about the gas temperature at the outer radius. This 
temperature is typically set to the Compton temperature assuming that the central X-rays 
heat a gas up to an equilibrium Compton temperature. However, in some of our simulations 
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the optical depth toward the radiation source or the local density in the flow is so high that 
radiation cannot heat a gas at large radii to the Compton temperature. In such cases, there 
is a mistmatch between the gas temperature assumed at To (i.e., Tq) and that computed for 
r close to Tq. 

It is possible that gas at large radii is heated not only by the central source but also by 
shocks or other sources such as supernovae. To mimic such sources of heating, we introduce 
an uniform background X-ray radiation, J-'x,h in some of our simulations. We illustrate 
effects of this radiation by showing results of one model (run Crbgd), where we assumed 
J-'x,h = 1-2 X 10'' erg cm~^ s~^, for which the gas with relatively low density is comptonizcd, 
e.g., for p = 1 X 10~^° g cm""^, ^ = 2.5 x 10"^ (we assume that Tr of the background radiation 
is the same as of the central source, that it 8 x 10^ K.) 

3. Results 

As in Paper I, we assume the mass of the nonrotating BH, Mbh = 10^ M0 and the 
disk inner radius, = 3 rs = 8.8 x 10^^ cm throughout this paper. We compute the 
total accretion luminosity as L = r/MaC^ = 2riGMBRMg,/rs, where r] and Ma is the rest mass 
conversion efficiency and the mass accretion on BH, respectively. We assume a relatively high 
conversion efficiency appropriate for disk accretion onto a non-rotating BH, i.e., r) — 0.0833. 

We express the accretion luminosity in in units of the Eddington luminosity for the 
Schwarzschild BH, i.e., Lsad = ^ncGMBn/cre- We refer to this normalized luminosity as the 
Eddington number, F = L/L-^dd = (<7e^a)/(87rcrs). 

Table 1 summarizes the properties of models from Paper I (runs A, B, Bl, B2, B3, and 
C) and our new models (the other models listed in the table). Columns (2) to (11) give the 
input parameters that we varied: the Eddington number, F, the disk contribution to the total 
luminosity, /disk, the corona contribution to the total luminosity, the UV contribution to 
the total luminosity, /uv, the X-ray contribution to the total luminosity, fx, the radiation 
temperature, Tr, the X-ray background radiation flux, J^x,h, the gas temperature at the outer 
boundary, Tq, the gas density at the outer boundary, pq, and the circularization radius, r^. 
Columns (12) to (17) give some of the gross properties of the solutions: the mass inflow 
rate through the outer boundary Min(ro), the net mass flux rate through the inner boundary 
-^net(^i); the mass outflow rate through the outer boundary Mout(^o); the maximum outflow 
velocity at the outer boundary, v-^, the outflow power carried out through the outer boundary 
in the form of kinetic energy, P]^{ro), and in the form of thermal energy, Pth('"o)- Table 1 
also explains our convention of labeling our runs. All other model parameters not listed in 
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Table 1 are as in Paper I. 

3.1. Effects of Gas Rotation 

Simulations without rotation presented in Paper I, show that an infaUing gas coUimates 
an outflowing gas and that the coUimation increases with increasing radius. In addition, for 
a given F, the coUimation increases with decreasing ratio between fx and /uv- Out of three 
cases explored in Paper I, case C with the smallest /uv//x, shows the strongest coUimation 
and the highest efficiency of turning an inflow into an outflow. Fig. 1 in Paper I and Fig. 1 
here show that in run C the gas is siphoned off within a very narrow channel along the pole. 

In Paper I, we argued that the coUimation, outflow power and other results will likely 
change if one would allow for signiflcant gas rotation. One would expect that the gas will 
converge toward the equator due to the combination of the centrifugal and gravitational 
forces. This, in turn, will likely broaden and weaken the outflow in the polar region because 
less gas will be pushed toward the polar region. 

To test effects of rotation, we rerun models presented in Paper I with rotation. We set 
— 300 that is near the maximum value for which the flow will not circularize inside our 
computational domain. We do not consider higher at the moment because we want to avoid 
complexities that will result from formation a rotationally supported torus or disk inside 
the computational domain (e.g., Hawley, Smarr & Wilson 1984a; 1984b; Clarke, Karpik 
& Henriksen 1985; Hawley 1986; Molteni et al. 1994; Ryu et al. 1995; Chen et al. 1997; 
Toropin et al. 1999; Kryukov et al. 2000; Igumenshchev & Narayan 2002; Proga & Begelamn 
2003; Chakrabarti et al. 2004, and references therein). In these exploratory simulations, our 
choice of high r^, yielding low I, allows us then to study flrst, relatively simple flows and set a 
stage for modeling more complex flows with high I. We assume that the circularized gas will 
accrete onto SMBH on a viscous time scale. However, as we do not model this part of the 
flow, we do not consider any details of an actual process/cs leading to transport of angular 
momentum. The transport is most likely due to magnetorotational instability (Balbus & 
Hawley 1991) but a contribution from the photon viscosity can be important in the case of 
high radiation luminosity. 

We note that our choice of low / can be relevant to real QSOs because feeding a SMBH in 
QSO with gas of high / can be problematic. Namely, for high /, an accretion disk would form 
at large radii and be self-gravitating (e.g., Paczynski 1978; Shlosman & Begelman 1987). 
Converting disk material to stars could then starve the SMBH and QSO would be quenched 
(e.g., Goodman 2003). After reviewing variety of possibilities, Goodman (2003) suggested 
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that QSO disks do not extend beyond a thousand rs so that they could be gravitationally 
stable. If so such the disks must be replenished with gas of small I as that we explore here. 

Figure 1 presents the results for runs C, CR, and Cr. For run CR, a step function 
describes the angular distribution of angular momentum on the outer boundary (we set 
to 45° in see eq. 4). For 45° < 9,we assume that the specific angular momentum at the outer 
boundary equals Iq, whereas for 6 < 45°, / = 0. For run Cr, a smooth function describes I 
(see eq. 5). The figure shows the instantaneous density and temperature distributions, and 
the poloidal velocity field of the models. In addition, it shows the Compton radius corrected 
for the effects of radiation pressure due to electrons (sec eq. [19] in Paper I) and the contours 
of the Mach number equal to one, (M = ^yv^ + v^/cs = 1, where Cg = 'yP/p is the sound 
speed). 

The detailed calculations confirm our general expectations: compared with the non- 
rotating case, the outflow is less coUimated and weaker in the rotating case. As Table 1 
shows, in the runs with rotation (runs CR and Cr), the outflow power, Mout(^"o), and Vr are 
lower than those in the run without rotation (run C; compare also flg. 4 in Paper 1 with 
Figs. 2 and 3 here). Another difference is that in runs CR and Cr, gas does not cool as much 
as in run C, especially at small radii (i.e., r' < 1 x 10^). Comparison between runs C, CR, and 
Cr shows that in run CR the solution is an intermediate one between runs C and Cr which 
is not too surprising, given the step function is an intermediate distribution between zero-/ 
and the one described by function fi. In particular, the outffow in run CR is less coUimated 
than in run C and more coUimated than in run Cr. One of new unexpected features that we 
found in run Cr is that the relatively cold outffow is fragmented and time-variable. 

Although the ffow in run Cr settles down into a time-averaged steady state, it is not as 
steady as in run C. An indication of this behavior can be found in Fig. 3 that shows three 
radial mass ffow rates as a function of radius: the net rate, M^ct{r), the inffow rate Afin(r), 
and the outffow rate Mout{f^) (see eqs. 22, 23, 24 in Paper I for formal deffnitions). For a 
perfect steady state, one expects Mnet('") = Mn(?') + A/out (? ) =const at all radii as in run 
C (see Fig. 4 in Paper I). However, in case Cr the above equation holds only at small radii, 
r' ^ 10^ K. We note that contrary to Mout: is a smooth function of radius. Thus the 
unsteadiness of in the ffow appears to be caused by the unsteady behavior of the outffow, 
especially the outffow at large radii where it can cool down. 

We relate the fragmentation and time- variability of the outffow to line force turning on 
abruptly when T decreases below ~ 5 x 10^ K and turning off when T increases again. Fig. 1 
shows that T decreases in the regions where the inffow sharply turns into the outffow. The 
density there increases and the gas radiatively cools. The turning on of line force leads to an 
enhanced acceleration of the outffow, but this alone is not sufficient to fragment the outffow. 
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For example, in runs C and CR the cold outflow is not fragmented and is quite steady. Thus 
there must be another factor/s that may contribute to fragmentation and time- variability. 
We note that in run C, the outflow is nearly radial hence its inner parts shield the outer parts 
form the central radiation. Consequently, the outflow can not be heated downstream by the 
central radiation. However, in run Cr the outflow is not radial and the flow can be heated up 
downstream because, as its density decreases during acceleration, it is irradiated by stronger 
unattenuated X-ray flux. There the outflow orientation with respect to the radiation flux 
appears to be one of the key factors causing fragmentation and time- variability of the outflow. 
This conclusion is supported by the fact that even in run Cr, the outflow is not fragmented 
at large radii where the outflow becomes almost radial, and clumps merge with each other. 

To show the variable solution in more details, in Fig. 4 we present a sequence of density 
maps of the inner part of the flow in run Cr at four different times. The left panel shows 
the flow at a time when a clump brakes from a high density fllament at z' ~ 1.5 x 10^ r*. 
Subsequent panels show how this and other clumps move outward and how they are stretched. 
Fig. 4 shows also formation a new clump (the second panel from the right). We find that 
clumps form usually at the same location (i.e., at r' 5 x 10^ r* and 2;' 1.5 x 10"^ r*) 
every 10^^ s or so which is of order of a dynamical time scale at radius where the clumps 
form. Generally, despite the time-variability, the instantaneous maps shown in Fig. 1 are 
representative of run Cr because they show an example of a large scale inflow and outflow 
with continuous production of small scale clumps that merge at large radii (i.e., beyond 
r' > 10^). 

Figs. 3 and 4 show that in runs CR and Cr the outflow power is dominated by the kinetic 
energy not the thermal energy. However, comparing with run C (see Fig. 4 in Paper I) the 
dominance is not as strong. 

We conclude that in case C, rotation reduces the outflow collimation and the outward 
flux of mass and kinetic energy. Rotation also leads to fragmentation and time- variability of 
the outflow and an increase of the outward flux of the thermal energy. As expected, rotation 
does not change much the mass inflow rate through the outer boundary. 

Figs. 5 and 6 show results for case B with and without rotation, i.e., runs B and Br 
(see also Fig. 3 in Paper 1). In this case, rotational effects are almost the same as in case C. 
The main difference is that in run Br, an outflow does fragment and the overall flows settles 
down into a steady state. This however is understandable: in run Br, radiative heating is 
strong and the gas does not cool therefore line force does not turn on. 



- 11 - 



3.2. Effects of the position-dependent radiation temperature 

We return now to case C and consider effects of the position-depenedent Tr. Left panels 
of Fig. 7 and Fig. 8 show resuhs for run Cx. In comparison with run C, the outflow in run 
Cx is broader. The mass outflow rate in run Cx is only slightly higher than in run C. 
However, the outflow rate almost cancels out the inflow rate so that the net rate is two 
orders of magnitude smaller than the mass flux through the outer boundary. Run Cx is a 
good example where AGN irradiation can signiflcantly reduce the rate at which the central 
engine is fed with matter. 

In run Cx, Tr is lower near the poles than near the equator. In addition, Tr near the 
pole is lower in run Cx than Tr in run C. The latter difference explains why the outflow in 
run Cx is so strong: the relatively low 7r in the polar region leads to a lower gas temperature 
in this region. This in turn leads to more mass being pushed towards the pole. This mass 
can then be effectively turned into an outflow because in the polar region the radiation flux 
is sub-Eddington and additionally, line-force turns on there because the gas temperature is 
low enough. In other words, the siphon effect seen in many of our simulations is very strong 
in run Cx. 

Right panels of Fig. 7 and Fig. 9 show results for run Crx that is a rerun of run Cx with 
rotation. Comparing these two run, we flnd that effects of rotation in the simulations with 
the position-dependent Tr are similar to those in the simulations with constant Tx. Namely, 
rotation decreases the degree of the outflow coUimation and decreases the outward flux of 
mass and kinetic energy. In addition, rotation leads to an increase of the outward flux of 
the thermal energy. In run Crx, the cold outflow is nearly radial and does not fragment as 
much as in run Cr. 

In this paper, we do not present results for case A with rotation because in this case 
the flow is dominated by thermal effects and rotation does not change much the solution. 
However, we present here results for case A with the position-dependent Tr (run Ax) because 
they show new effects. 

Fig. 10 compares results for runs A and Ax. In run Ax, the flow pattern is different that 
that seen in runs presented in Paper I or runs shown here so far. The dramatic difference 
is that in run Ax the equatorial inflow is replaced by an equatorial outflow. Generally, 
in run Ax there are equatorial and polar outflows both being fed by an inflow of gas at 
intermediate polar angles. The equatorial outflow is a simple consequence of the higher 
radiation temperature near the equator that leads to a high gas temperature and an enhanced 
thermal expansion. In the polar region, where the gas temperature is lower, an outflow is 
driven by thermal expansion assisted by radiation pressure on electrons. 
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Overall the flow pattern in run Ax is dominated by the outflow. This leads to the net 
mass inflow rate at small radii being one order of magnitude lower than the inflow rate at 
large radii. 

We conclude that both rotation and position-dependent Tr lead to qualitative and quali- 
tative changes in the flow. Most prominent of them are weaker coUimation and fragmentation 
of the outflow in cases with rotations, and production of an thermal equatorial outflow in 
cases with position-dependent Tr. 

3.3. Complex Case 

A reahstic model of accretion flows should include many physical effects. Here, we 
focused on the role of gas rotation and position-dependent Tr. We consider our simulations 
just as exploratory tests. These tests support the notion that AGN radiation can play a very 
important role in determining the ionization structure and dynamics of matter on sub- and 
parsec-scales. We finish the presentation of our models with showing results for two runs 
that illustrate how complex the flow dynamics can be even in a very simple set-up as the 
one we focus on here. 

One of our motivations is to understand gas dynamics in broad line regions (BLRs) and 
Narrow Line Regions (NRL) so characteristic to AGN. These regions are thought be made 
of cold gas clouds moving randomly or nearly randomly and having a small flUing factor 
(e.g., Krolik 1999 and reference therein). Formation, evolution and other key aspects of 
these clouds are not well understood. Our simulations show that an accretion flow which is 
initially smooth and spherical can break into inflows and outflows. In cases with rotation, 
we have seen outflows fragmented due to hue force and X-ray irradiation. These results raise 
the following question: can we produce many cold clouds with a small flUing factor. The 
answer seems to be yes as we show in Figs. 12 and 13, and 14. 

Figs. 12, 13, and 14 compare results for case G with F = 0.9 and Po — 10~^° g cm"^. 
Left panels in Fig. 12 show results for run Grgd. The density and temperature distribution 
and as well as some other properties of this run differ a lot from other runs shown here. The 
main difference is a much larger dynamical range in the temperature and density plus the 
fact that the ffow is far from reaching any steady state. However, close inspection of the 
results from run Grgd show that the ffow pattern is similar to that seen in most other runs: 
there is an equatorial inffow and a bipolar outffow. The latter is not fragmented because it is 
radial in agreement with our explanation of the outflow fragmentation. The large dynamical 
range in the density and temperature in run Grgd is a simple consequence of the higher 
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density at the outer boundary. A high density gas can cold much faster than a low density 
gas. In this run the gas is also heated by shocks because for F = 0.9 radiation pressure of 
electrons alone can drive a powerful and broad outflow that collides with an inflowing gas. 

As we discussed in §2.2, an accretion flow at large radii does not have to be heated 
by the central radiation source. The lower left panel in Fig. 12 illustrates this point as one 
can sec that near the equator a good fraction of the inflow is relatively cold. One of the 
prediction of such a model is that for a wide range of the inclination angle, this cold gas 
should produce absorption lines redshifted with respect to the systematic velocity. However, 
such lines are not being observed. 

The problem of large scale accretion of cold gas can be easly overcome by introduction 
of non-central source of heating. Right panels in Fig. 12, show results for run Crbgd which 
is a rerun of model Crgd with the X-ray background radiation (see §2.2). 

In run Crbgd, there are no large regions of shock heated gas. The only region where 
shock heating is important is in a narrow polar region of a low density. The background 
radiation heating in this run keeps the gas from rapid cooling which in turn can lead to 
abrupt turn-on of line driving and strong expulsion of gas as seen in run Crgd. 

Comparing runs Cr and Crbgd, we sec that an outflow tends to fragment more if the 
gas density is higher. The density and temperature maps show that the outflow occupies 
a relatively large fraction of the computational domain. However, Fig. 14 shows that this 
outflow docs not change much the overall mass budget. As in run Cr, this is caused by 
rotation that reduces the amount of gas that is pushed toward the polar region where it can 
be siphoned off. 

4. Conclusions 

We have calculated a series of models for rotating flows that are under the influence of 
supermassivc BH gravity and radiation from an accretion disk surrounding the BH. We seek 
to determine sclf-consistently what fraction of the flow is gravitationally captured by the BH 
or what fraction is driven away by thermal expansion and radiation pressure. This work is 
an extension of the work presented in Paper I, where nonrotating flows were studied. Here, 
we consider effects of rotation and of a position-dependent radiation temperature, density 
at large radii, and an uniform X-ray background radiation. 

As in the non-rotating case, the rotating flow settles into a conflguration with two 
components (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving 
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the system along the pole. However, the rotating flow does not reach a steady state. In 
addition, rotation reduces the outflow coUimation and the outward flux of mass and kinetic 
energy. Moreover rotation increases the outward flux of the thermal energy and can lead 
to fragmentation and time-variability of the outflow. In future, we plan to check whether 
thermal instability can contribute to fragmentation and time-variability of the outflow. As 
expected, rotation does not change much the mass inflow rate through the outer boundary. 

In our model, the radiation comes from an UV emitting disk and X-ray emitting spherical 
corona. As a result, the radiaton temperature is position-dependent: in the polar region 
radiation is dominated by a softer disk component whereas near the equator, radiation 
is dominated by a harder corona component. The two main changes due this position- 
dependence are: (i) an increase in the power of the outflow in the polar region and (ii) 
development of a large scale thermally driven outflow in the equatorial region. 

Overall, we conclude that our exploratory study provides an additional support to the 
idea that AGN radiation can signiflcantly change gas dynamics and photoionzation structute 
on sub-parsec- and parsec-scales. As it have been discussed and shown in the past, we found 
that AGN radiation can significantly reduce the rate at which the central BH is fed with 
matter (e.g.. Figs. 8 and 13). This result should not depend on the inner radius of the 
computional domain because most of the outflow is launched from a radius larger than the 
inner radius of the computational domain. However, we note when reducing the inner radius 
of the computional domain one should also consider additional processes, in particular, disk 
accretion and disk winds and jets. Thus our mass inflow rate should be viewed only as an 
upper limit for the BH accretion rate. 

Our results emphasize also a little appreciated feature, i.e., AGN radiation can drive 
a non-spherical, multi-temperature and very dynamic flow pattern. This result may have 
implications for the problem of AGN feedback and the problem of the origin, geometry and 
physics of NLR and BLR. 

This work is supported by NASA through grants HST-AR-10680 and HST-AR-11276 
from the Space Telescope Science Institute, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Table 1: Summary of results 
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Crx 
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0.05 
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-11 
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10 


0.01 


Crgd 
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0.05 
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Note. — (1) in 2 X 10^ K, (2) in 10^ org cm-^ s-\ (3) in lO^^i g cm--\ (4) in (5) in lO^-'^ g s-\ (6) 
in km s~^, (7) in 10'*° erg s~*. We use the following convention to label our runs: the first character refers 
to the values of /x and /uv: A is for /uv = 0.5 and /x = 0.5, B is for /uv = 0.8 and /x = 0.2, and C is 
for /uv = 0.95 and /x = 0.05. Runs A, B, and C, are the fiducial runs. If the first character is followed by 
a lower case letter or letters or a number, it means that it is the same run, but modified by introduction 
of: rotation (letter 'R' and 'r' if we use eq. 4 or 5, respectively), the position-dependent Tr (letter 'x'), the 
X-ray background radiation (letter 'b'), a higher F (letter 'g'), and a higher po (letter 'd'). The numbers 1, 
2, and 3 correspond respectively to Tr = 1/10, 1/3, 3 in units of 2 x lO'^ K. 
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Fig. 1. — Comparison of the results for runs C, CR, and Cr (left, middle, and right column, 
respectively) Top row: Maps of logarithmic density overplotted by the poloidal velocity. 
For clarity, the arrows are plotted with the maximum velocity set to 1000 km s~^. The 
solid curves show are the contours of the Mach number equal one. Bottom row: Maps of 
logarithmic temperature overplotted by the direction of the poloidal velocity. The solid curve 
in the bottom left corner marks the position of the Compton radius corrected for the effects 
of radiation pressure due to electron scattering (see eq. 19 in Paper I). The length scale is 
in units of the inner disk radius (i.e., r' = r/r^ and z' = z/r^). 
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Fig. 2. — Top panel: The mass flux rates as a function of radius for run CR. The sohd, 
dashed, and dotted hne corresponds to the outflow, inflow, and net rates, respectively (see 
eqs. 20, 21, and 22 in Paper 1 for the formal deflnitions). Note that the absolution value 
of the inflow and net rates are plotted because these quantities are negative. Bottom panel: 
The energy fluxes carried out by the outflow as a function of radius in run CR. The solid 
and dashed line corresponds to the kinetic and thermal energy flux, respectively (see eqs. 23 
and 24 for the formal deflnitions). The length scale is in units of the inner disk radius (i.e., 
r' = r/r*). 



-20- 




- 21 - 




Fig. 4. — A sequence of density maps for the inner part of run CR after 4.95, 5.00, 5.05, and 
5.10x10^2 s (from left to right). Run CR is an example of unsteady flow discussed in detail 
in section 3.1. As in fig. 1, the length scale is in units of the inner disk radius (i.e., r' = r /r^, 
and z' = z/r^). However, the r' and z' ranges are 2.5 time smaller compared to fig. 1. 
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Fig. 5. — As in Fig. 1, but for runs B and Br. 
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Fig. 6. — As in Fig. 2, but for run Br. 
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Fig. 7. — As in Fig. 1, but for runs Cx and Crx. 
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Fig. 8. — As in Fig. 2, but for run Cx. 
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Fig. 10. — As in Fig. 1, but for runs A and Ax. 
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Fig. 12. — As in Fig. 1, but for runs Crgd and Crbgd. 
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Fig. 13. — As in Fig. 2, but for run Crgd. Note an increase range along the y-axis compared 
to Fig 2. 
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Fig. 14. — As in Fig. 2, but for run Crbgd. Note an increase range along the y-axis compared 
to Fig 2. 



